Search for a T-odd, P-even Triple Correlation in Neutron Decay 
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Background: Time-reversal-invariance violation, or equivalently CP violation, may explain the 
observed cosmological baryon asymmetry as well as signal physics beyond the Standard Model. In 
the decay of polarized neutrons, the triple correlation D{Jn) ■ (Pe xpv) is a parity-even, time-reversal- 
odd observable that is uniquely sensitive to the relative phase of the axial-vector amplitude with 
respect to the vector amplitude. The triple correlation is also sensitive to possible contributions from 
scalar and tensor amplitudes. Final-state effects also contribute to D at the level of 10~^ and can be 
calculated with a precision of 1% or better. Purpose: We have improved the sensitivity to T-odd, 
P-even interactions in nuclear beta decay. Methods: We measured proton-electron coincidences 
from decays of longitudinally polarized neutrons with a highly symmetric detector array designed 
to cancel the time-reversal-even, parity-odd Standard-Model contributions to polarized neutron 
decay. Over 300 million proton-electron coincidence events were used to extract D and study 
systematic effects in a bhnd analysis. Results: We find D = [— 0.94± 1.89(stat) ±0.97(sys)] x 10"*. 
Conclusions: This is the most sensitive measurement of D in nuclear beta decay. Our result can 
be interpreted as a measurement of the phase of the ratio of the axial-vector and vector coupling 
constants {Ca/Cv = |Aje'"*^^' ) with 4>av = 180.012° ±0.028° (68% confidence level) or to constrain 
time-reversal violating scalar and tensor interactions that arise in certain extensions to the Standard 
Model such as leptoquarks. This paper presents details of the experiment, analysis, and systematic- 
error corrections. 



PACS numbers: 24.80.-fy, ll.30.Er, 12.15.Ji, 13.30.Ce 
I. INTRODUCTION 



The symmetries of physical processes under the trans- 
formations of charge conjugation (C), parity (P), and 
time reversal (T) have played a central role in the de- 
velopment of the Standard Model of elementary-particle 
interactions [T]. Time-reversal-symmetry violation (or T 
violation), which is equivalent to CP violation assum- 
ing CPT symmetry, has been of particular interest be- 
cause it is sensitive to many kinds of new physics. The 
CP-violating parameters of the Standard Model are the 
Cabibbo-Kobayashi-Maskawa (CKM) phase, which en- 
ters in the mixing of three generations of quarks, and 
the parameter 6qcd- The effect of the CKM phase is 
strongly suppressed in the permanent electric dipole mo- 
ments (EDMs) of the neutron and heavy atoms O Hj , 
and recent EDM results combine to set upper limits on 
(^QCD- AH laboratory measurements to date are consis- 
tent with a single source of CP violation, i.e. the phase 
in the CKM matrix. An exception may be the 3.2 sigma 
deviation observed recently as an asymmetry in the pro- 
duction of pairs of like-sign muons reported by the DO 



collaboration [S]. 

In spite of this success, laboratory and astrophysi- 
cal observations, which include neutrinos with non-zero 
masses, the abundance of non-baryonic dark matter, and 
the baryon asymmetry of the universe, provide strong 
evidence that the Standard Model is incomplete. Gen- 
eration of the baryon asymmetry requires CP viola- 
tion that cannot be accounted for by Standard-Model 
physics [HI [7] . This provides strong motivation to search 
for new sources of CP violation. Such CP-violating 
phases would also, in general, affect T-odd observables 
in neutron decay, in particular, the T-odd/P-even triple 
correlation in polarized neutron decay. 

We have measured the triple-correlation in the decay 
of polarized neutrons at the NIST Center for Neutron 
Research (NCNR) using the emiT apparatus [SHTO]. The 
D-coefficient is uniquely sensitive to the relative phase of 
vector and axial vector amplitudes and is also sensitive 
to scalar and tensor currents. Our result, first reported 
in reference fll , sets a new limit on this phase as well as 
on certain combinations of scalar and tensor currents. 

This paper presents the details of the experiment and 
analysis and examines the implications of this result. The 
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paper is organized as follows: The context of the mea- 
surement is presented in the remainder of this introduc- 
tion. In section [IT] the measurement principle based on 
the symmetries of the apparatus is presented. The Monte 
Carlo simulations used to understand the apparatus and 
important details of the performance of the apparatus are 
provided in sections |III| and |IV| In section [V] we present 
the details of the data set; and section [Vl| is a detailed dis- 
cussion of all systematic effects and corrections applied 
to the data. In section [VII| we describe the principle of 
the blind analysis and present the final result; and in 
section |VIII| we summarize the implications of our result 
and prospects for improved measurements. 



A. Polarized neutron decay 

Neglecting recoil-order corrections, the general form of 
the Hamiltonian for beta-decay {e.g. n — >■ peP) can be 
written [T^ 



Hint — 



1 
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h.c 



(1) 



where the subscripts V , A, S, and T respectively re- 
fer to vector, axial-vector, scalar, and tensor contribu- 
tions, and we have left out the pseudoscalar amplitude, 
which vanishes in the limit of non-relativistic nucleons. 
Allowing the possibility of T violation, the C's are com- 
plex numbers, and there are 19 free parameters plus 
a single arbitrary overall phase. The Standard Model 
is written with only left-handed V and A interactions 
with Cv — C'y and Ca — C^. Thus the number of 
free parameters is reduced to three: |Cy|, \Ca\ and 
the relative phase of Ca/Cv = lAje*"*^^. The value 
of \Cv\ — Gp\Vud\ follows from the conserved- vector- 
current hypothesis (CVC) with Gp determined from the 
muon lifetime, and \Vud\ most precisely determined from 
super-allowed beta decays [13]. The parameter |A| is de- 
termined from other measurements including the beta 
asymmetry (A-term) and electron-neutrino correlation 
(a-term) . 

The Hamiltonian of Eq.[l] leads to the following differ- 
ential decay rate for polarized neutrons and no final-state 
electron polarization [T^]: 
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where S{Ee) = F{Ee)peEeiE""'=' - E^f is the phase 
space factor with F{E^) the Fermi-function for Z = 1, 
Pe(£'e) and TpviEv) are the momentum(energy) of the 



electron and antineutrino, and the neutron polarization 
P = -^y^ is the ensemble average of the neutron spin. 
The triple-correlation, I?(J„) • (pe x Pi,) is P-even but 
odd under motion reversal. Time-reversal is the com- 
bination of motion-reversal and initial/final-state rever- 
sal. Thus contributions to D can originate from T- 
violating interactions and from final state effects, i.e. 
D = D'^+Dpsi [m, where 



Dfsi ~ 1-1 X 10" 



Pe 



0.3 X 10 



n" 

-5 Pe 



(3) 



For the neutron Dpsi ~ 1-2 x 10 ^ and can be calculated 
to 1% or better [TMI]- 

In the neutron rest frame, P;, = — (Pp + Pe) (for cold 
neutrons, the neutron velocity in the lab frame has a 
negligible effect), and the triple correlation can be written 



D{3n) ■ (pe X p^) = D{in) ■ (Pp X Pe). 



(4) 



We therefore extract D by measuring proton-electron an- 
gular correlations in polarized-neutron decay. 



B. The Physics of 

In the context of Eq. [T | £) y depends on the values of 
Ci and can be written as jl5j 

= -^\Mf\\Mgt\IHCvC*a + C'yC'X] 
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where \Mp\ and |Mgt| are matrix elements for 
Fermi and Gamow-Teller transitions (for neutrons, 

\Mgt?/\Mp\'' = 3), and 

^ = |M^P(|Cvp + |Ct.p + |C5p + |C^P) 

+ \Mgt?{\Ca? + \C'a? + \Ct? + \C't\^). (6) 

We assume terms quadratic in Cs and Ct can be ne- 
glected in ^, and take Cy — C'y and Ca = (i-e. 
left-handed V-A interactions) so that 
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where A = lAle*"^-*^ = Ca/Cv- The first term is sensi- 
tive to the phase of A and can be approximately written 
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D^^ « 0.435 sin .^AV- for |A| = 1.2694 [H]. The remain- 
ing terms show sensitivity to combinations of scalar and 
tensor amplitudes that signal beyond-Standard-Model 
Physics. 

Standard-Model CP violation makes very small con- 
tributions to Df. The CKM phase enters light-quark 
processes at loop level yielding sa IQ-^^ ggjj^g The 
parameter Oqcd is tightly constrained by the EDMs of 
the neutron and ^^^Hg and also makes a very small con- 



tribution to Dn 



I.e. 



less than 10"" 



Beyond- 



Standard-Model physics, in particular left-right symmet- 
ric models, exotic fermions, and leptoquarks give rise to 
couplings that contribute to at a level comparable to 
the current experimental sensitivity; however, it has been 
argued that limits on T-odd/P-odd observables, in par- 
ticular EDMs, can be used to place limits on T-odd/P- 
even interactions and thus on Dj by calculating the T- 
odd/P-even effects of radiative loops on T-odd/P-odd 
interactions |18j-t22j. For the neutron EDM, this limit 
is more stringent than the sensitivity of the experiment 
described here by as much as an order of magnitude; 
however, this argument is based on the assumption of 
complete absence of cancellations between different con- 
tributions to the neutron EDM, which cannot be a priori 
excluded [22 . Table |l] summarizes the contributions to 
Dj from the Standard Model and extensions. 



TABLE L Expected contributions to Z)^ for neutron de- 
cay from parameters of the Standard Model and beyond- 
Standard-Model physics based on measurements in other sys- 
tems. The broad range of limits arises in the cases of signifi- 
cant model dependence. 



Model/Parameter 



CKM phase 
Oqcd 

Left-right symmetry 
Non-SM fermions 
Charged Higgs SUSY 
Leptoquark 
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Limit on 

2 X 10 
10-^-10" 
lO-'^-lO" 
10-^-10" 
10^^-10" 



C. Recent results 

The two most recent measurements of D 
in neutron decay are the emiT-I measurement 
and the TRINE experiment result from Insti- 
tute Laue Langevin, Grenoble (ILL). For emiT-I, 
D=[-%± 12(stat) ± 5(sys)] x 10"" and the TRINE 
result was D ^ [-2.8 ± 6.4(stat) ± 3.0(sys)] x 10"'' [53]. 
The current Particle Data Group average for the neu- 
tron, D — [—4 ± 6] X 10""' PU, also includes results 
from references p5H?7] . A measurement in ^^Ne, where 
the final-state-interactions are more than an order 
of magnitude larger than for the neutron, resulted in 
Di9No = [0.7±6] X 10-4 inOlISI]- We also note that the R 
coefficient of the T-odd/P-odd correlation J„ • (pe x cTg), 



which is linearly sensitive to S and T amplitudes, has 
recently been measured for the neutron |32j and for 

^Li m. 



II. THE EMIT-II EXPERIMENT 

The emiT experiment was designed to measure proton- 
electron coincidences in the decay of neutrons polarized 
along the axis of an array of detectors. The symmetry of 
the detector array allowed us to discriminate the triple 
correlation from the T-even/P-odd A- and B-coefScient 
correlations. The layout of the experiment is shown in 
Fig. [T] The cold neutron beam was transported by the 
neutron guide NG6 to the experiment. The neutrons 
were polarized and passed through a spin flipper. Down- 
stream of the spin-flipper, neutron spins were adiabati- 
cally transported through rotation of the magnetic field 
to longitudinal, along the axis of the detector array. 

The detector array, illustrated in Fig. [2] consists of four 
electron-detectors alternating with four proton-detector 
planes arranged in an octagonal geometry concentric 
with the neutron beam. Each of the four proton-detector 
planes consists of 16 separate cells arranged in two rows 
of eight cells. The protons are detected by negatively- 
biased surface barrier detectors (SBDs) that are incor- 
porated into a focusing cell illustrated in Fig. [3] Within 
the fiducial volume of the detector array, neutrons are 
polarized parallel or antiparallel to the magnetic field de- 
pending on the state of the spin flipper. The magnetic 
field in the detector region is approximately 560 fiT in 
magnitude and is aligned parallel to the neutron beam 
and detector axis. 



A. Electron-proton Coincidence Events 

The data set consists of 512 sets of coincidence events 
from the combination of 64 proton cells and four electron 
detectors for the two spin-flipper states. The total num- 
ber of counts for a given run time is labeled as N^'^\ 
where the ± indicates the spin-flipper state (neutrons 
nominally parallel/opposite to the magnetic field B), pi 
labels the proton-cell, and tj labels the electron detector. 
The neutron-spin dependence of the count rates depends 
on the correlations given by A, _B, and D and is given 
by the difference of rates for the two spin-flipper states, 
while the average of rates includes the spin-independent 
beta-neutrino correlation (a term). 

In order to isolate the neutron-spin dependent terms, 
we deflne the asymmetry 



(8) 



In principle, the N^'^' follow from integrating Eq. [2]over 
the neutron beam, the detectors' acceptances, electron 
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FIG. 1: Layout of the emiT-II experiment beam-line. The neutron beam is nominally unpolarized upstream of the polarizer, 
and is vertically polarized downstream. Polarized neutrons are guided within beryllium coated glass tube to the detector. As 
shown, the spin flipper reverses the direction of the vertical magnetic field B over a short distance so that the neutron spin 
Jn, which remains polarized vertically upward, reverses with respect to the magnetic field. Downstream of the spin flipper, 
solenoids rotate the magnetic field into the horizontal direction, parallel to the neutron beam. 



momentum, and neutrino angles for a fixed time so that where 0^ and 0^ are the polar angles of the electron and 

antineutrino with respect to the neutron polarization P, 
and 



A(/?eP-Pe)+i?(P-p.)+i?(/?e(f^)P-(Pp X p^)) 



(1) + a(/3ePe ■ Pi/) 



(9) 



where /3e — We/c is the electron velocity, P is the neutron 
polarization at a given position, and the brackets (( )) 
indicate that each term is integrated over energies, the 
neutron-beam distribution, and sohd angles for proton 
detector pi and electron detector ej . As shown below, the 
Z3-coefficient term can be isolated by forming a specific 
combination of w^'^j that cancels the parity- violating A 
and B correlations. 



B. The ideal experiment 

In order to understand the analysis technique, we begin 
by considering an ideal experiment with uniform longitu- 
dinal neutron polarization (P = -Pz), uniform neutron- 
beam density, and uniform efficiencies for all proton and 
electron detectors. Consider a proton detected in pi in 
coincidence with an electron detected in either detector 
62 or 63 as shown in Fig. [2j (Coincidence rates of pi 
with 6i and 64 were about 15-25 times lower and were 
not used to extract D in this analysis.) Thus, for longi- 
tudinal neutron polarization, the asymmetry from Eq. [9] 
can be written 



w 



P.- 



+A{l3e cos( 



(Pp X Pe)) 

B(cos6,,)], 



(10) 



(1) + a(/3ePe • Pt,) 

It is useful to define an instrumental constant that char- 
acterizes the sensitivity of w''^'^^ to the triple correlation, 
i.e. K^^"' = dmP^^'/dD or 



«P'^^(/3e(^)z-(PpXpe)). 



(11) 



The K'^'^^ used in the analysis were determined by Monte 
Carlo studies and are discussed in section IVI Fl 

In order to isolate the triple correlation, we note 
that the longitudinal component of the cross product, 
z • {pp X Pe), has opposite sign for ^163 and P162 coinci- 
dences: for J5163, it is positive, and for pie2, it is nega- 
tive. Thus we form a difference of spin- flip asymmetries 
for the two electron detectors. For exa mple , for pi we 



have: v^^ 



■wP^"^). Using Eq. 



longitudinal polarization are 

+ P^[Kfi^-M/3eC0S( 

+ p|k^^Mcos0.) 
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the yP' for 



\p1e3 



{f3e COS t 



\P1e2] 



(cos 9,, 



(12) 



where Kd = U^d"'' ~ -^d"")- Due to strong anti- 
correlation of proton and electron momenta, the asym- 
metries depend strongly on the axial position of the pro- 
ton cell. (Data are shown in Figs. 12 and [13]) For ex- 
ample, assuming D = 0, the vP^ (and wP^^j are equal 
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FIG. 2: The emiT-detector array. Left: side view showing proton-detector planes, with 16 proton cells (2x8) in each plane, 
and 50 cm long electron detectors. Right: end view showing the four proton-detector planes and four electron detectors. The 
magnetic field, directed parallel to the average neutron velocity, causes the proton and electron trajectories to be curved as 
indicated by the greatly exaggerated paths shown. In this paper, we refer to the magnitude of the relative angles between a 
proton detector and electron detector, e.g. for proton detector pi, these are 35°, 55°, 125°, and 145°, for a, ei, 63, and €2 
respectively. For reference, distances are measured from the center of the detector array. 




FIG. 3: Left: proton-cell layout with calculated proton paths 
showing focussing properties of electrode geometry. Right: 
Monte Carlo generated hit pattern for neutron spin parallel 
(red) and opposite (blue) z. 



but opposite for the upstream proton cell and an axially 
symmetric downstream proton cell for a uniform neu- 
tron beam, e.g. v^^ — —v^^^. Also note that in the 
absence of the 560 /iT magnetic field, the w^' are op- 
posite for the adjacent proton cell, e.g. v^^ = —v^^. 
Thus, for an ideal experiment with uniform longitudi- 
nal polarization and beam, the average of u^' from an 
upstream-downstream pair of proton cells {e.g. v^^ and 
wPi^) or adjacent cells {e.g. v^^ and v^^) will cancel the 
beta-asymmetry and neutrino-asymmetry terms leaving 
the £)-coefhcient term; however, the magnetic field affects 
the average Pe differently for odd {e.g. pi) and even {e.g. 
P2) proton detectors. Thus we combine data for all four 
proton cells in order to isolate D. 



Transverse neutron polarization 

A small misalignment of the magnetic field with re- 
spect to the detector axis gives rise to transverse- 
polarization effects. In the limit of a small beam diam- 
eter and a small angular acceptance of the proton cell, 
the average proton momentum (pp) is the same for the 
two electron detectors {e.g. 63 and 62), and for transverse 
polarization, indicated by the subscript T, we can write 

« DTsin((^P - (j)P^), (13) 

where (pp in the azimuthal angle of the polarization, (fP^ 
is the effective azimuthal position of the proton cell as 
indicated in Fig. [2] and vt ~ O.46Fsin0p. The factor 
0.46 is consistent with both Monte-Carlo simulations and 
the transverse-polarization calibration runs discussed in 
section [VI E[ Thus, for uniform neutron density and po- 
larization, the sinusoidal dependence of on the az- 
imuthal position of the proton-cell averages to zero when 
data from all four proton-detector planes are combined. 



Extracting D 

In order to cancel the effects of transverse polarization, 
the beta and neutrino asymmetries, and the 560 ^T mag- 
netic field, we must combine the v^^ from at least 16 pro- 
ton cells symmetrically located with respect to the center 
of the detector: two adjacent cells upstream of the de- 
tector center paired with two adjacent downstream cells 
from all four proton planes. Each set of 16 proton cells 
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have the same l^;''*], i.e. z^' = ±2 , ±6, ±10, and ±14 cm. 
For example, the shaded detectors in Fig. [5] correspond 
to z^*' = ±10 cm. For each \zP^\ we define a measured 
quantity D given by 



D = 



1 



(14) 



\zPi |— const 



where — 0.378 ± 0.019 is discussed in section VI F 
Due to the symmetry with respect to the detector center, 
each of the four possible D is an independent measure- 
ment of the same nominal quantity. This provides cross 
checks and maximum statistical power. The D are sub- 
ject to a variety of systematic effects that are estimated 
and applied as corrections in order to dete rmine D. Cor- 
rections to D are discussed in section VI 



to evaluate the response, energy loss, and build up of 
condensation on the SBDs. 

The propagation of neutron spin components from the 
polarizer and through the spin flipper and guide field 
were modeled by numerical integration of the Bloch equa- 
tions using Monte Carlo techniques to generate neutron 
trajectories 



IV. APPARATUS PEFORMANCE 

Detailed descriptions of the experiment components, 
including beam line, polarization, spin flipper, spin- 
transport magnetic fields and the detector array for the 
emiT-II run are described in detail in references [HI fTOl 
l55tilT] . Here we summarize specific features that impact 
the analysis and systematic effects described in this pa- 
per. 



III. MONTE CARLO SIMULATIONS 



Monte Carlo simulations were used both in the design 
of the experiment and to estimate several of the system- 
atic effects that impact the analysis. Separate simula- 
tions addressed the detector responses, proton focusing 
cells, and neutron spin transport. 

To track electrons, we used Penelope [34 , which has 
been experimentally validated for a variety of kinematics 
of relevance to neutron decay [35] . To track protons we 
used custom code embedded into Penelope. The Monte 
Carlo code used the measured beam distributions and 
took into account the 560 fiT axial magnetic field and 
details of the detector as described below. 

Simulations used a detector model that included all 
surfaces visible to decaying particles including cryo- 
panels and the proton-cell ground planes. For simplicity 
in backscattering studies, the grids on the proton de- 
tectors were replaced in the model with solid foil, and 
the simulations were verified for this substitution. The 
PENELOPE-based Monte Carlo code was used to analyze 
the effects of the electron response function and backscat- 
tering of protons and electrons. Section VI B provides 



details on how the backscattering events were identified. 
Section |VI C| gives details on the issues associated with 
the scintillator response. 

A separate simulation using SIMION was developed 
to track proton trajectories within the proton cells for 
proton events generated by the PENELOPE-based Monte 
Carlo. The detailed proton-cell geometry was used to 
determine the electric fields and to analyze focusing ef- 
ficiency as a function of incident proton momentum and 
position on the focusing-cell grid. In order to minimize 
computer time, protons that penetrated beyond the grids 
were assumed to be detected with efficiencies calculated 



with these separate proton-cell simulations. Section |VIB 
describes issues related to proton backscattering. 

The energy-loss code SRIM 37\ was used to provide 
proton energy loss backscattering probabilities. In ad- 
dition, SBD models were used with SRIM calculations 



A. Beam 

Lead-backed lithium-glass collimators defined the 
beam and reduced backgrounds at the detector. The 
beam dump was located two meters downstream of the 
detector and was also composed of lithium glass. A 3 mm 
hole in the glass allowed neutrons to pass through and 
the beam rate to be monitored continuously throughout 
the run. The mean neutron capture flux (the average of 
the neutron fiux weighted by the inverse neutron velocity, 
1/vn) was monitored with a fission monitor |43| . 

Estimates of the distribution of neutron density were 
used in the Monte-Carlo modeling of the experimen- 
tal and systematic effects. The cold neutron beam 
from NG6 was mapped by activating a 25 /im thick, 
natural-abundance dysprosium foil. Neutron absorption 
on ^^^Dy (28.2% abundance) follows the l/w„ law, and 
the foil is thin enough to provide a sufficiently accurate 
measure of the neutron density. The activated foil was 
subsequently laid on a beta-sensitive film, and the ex- 
posure was measured with an image reader with pixel- 
resolution of 200 /xm by 200 /xm [32] • Maps 18 cm up- 
stream and 18 cm downstream of the center of the detec- 
tor array are shown in Fig. [4] The dominant features of 
the maps are the expansion of the beam (by about 5 mm 
in radius over the 80 cm length of the detector array) 
and the horizontal shift of the beam, which is due to the 
properties of the polarizer. The beam expansion affects 
the longitudinal cancelations of the u^' from upstream- 
downstream combinations and leads to a systematic ef- 
fect discussed in section IVIDI 



B. Polarization and spin flipper 

The neutron polarizer is a supermirror bender 
(PSM) 44 . Neutrons are polarized parallel to the ver- 
tical magnetic field of the PSM. The spin flipper con- 
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FIG. 4: Neutron-beam distributions for positions 18 cm up- 
stream of the center of the detector (top) and 18 cm down- 
stream of the detector center (bottom). Contours show the 
film exposure relative to the maximum at the center of the 
upstream map. 

sists of two closely-spaced current sheets with horizontal 
currents. The magnetic field from the upstream current 
sheet was parallel to the PSM field, and the current in 
the downstream current sheet could be set parallel or 
antiparallel to the upstream current sheet. For antipar- 
allel current, the vertical magnetic field reverses direction 
over a distance of about 1 mm corresponding to a time 
of about 2 fis in the rest frame of a cold neutron. This 
time is short compared to the inverse Larmor frequency, 
and the neutron spin does not follow the magnetic field 
adiabatically. The result is that the neutron spin re- 
mains oriented in the original upward direction while the 
magnetic field reverses from up to down, thus reversing 
the projection of the spin with respect to the magnetic 



field direction. Small transverse neutron-spin compo- 
nents arise due to imperfect transitions. These transverse 
components precess around the magnetic guide-field as 
the neutrons move through the apparatus. Thus the po- 
larization direction of a decaying neutron depends on the 
vertex of the decay as well as the neutron's velocity and 
the magnetic field. The distribution of neutron velocities 
leads to azimuthal averaging of the transverse polariza- 
tion components. The 560 field was chosen so that 
the orientation of the neutron polarization P varied by 
less than 2 x lO^'^ across the detector. Downstream of 
the spin flipper, the guide field adiabatically rotates the 
neutron spin into the longitudinal direction and remains 
longitudinal throughout the detector array. 

The polarization was measured and mapped using a 
second polarizer as an analyzer downstream from the 
spin flipper |10j . A map of the estimated neutron po- 
larization, assuming a perfect spin flipper and identical 
polarization and analyzing power (P = Ap), is shown 
in Fig. [5] Though neither of these assumptions is ac- 
curate, this provides a lower limit on the neutron po- 
larization |10) . The increase in polarization from left to 
right is a property of the polarizer and leads to effects 
that are discussed in section [VI E[ When averaged across 
the beam, the lower limit on the neutron polarization is 
P > 0.91 (90% c.L), and we use P = 0.95 ± 0.05 in the 
analysis of D. The polarization and spin-flipper charac- 
teristics appeared to be quite stable over long periods, 
and the measured polarization was consistent with that 
measured for the first emiT run, which set the lower limit 
of P > 93% [9]. 



C. Magnetic fields 

Within the detector region, the magnetic field is main- 
tained by a set of eight 0.95 m diameter coils uniformly 
spaced over 2 m. The outermost coils had independently 
adjustable currents, while the inner six coils were con- 
nected in series. The current, alignment, and position 
of each coil were set to optimize the longitudinal field 
uniformity. An array of longitudinal coils was deployed 
as a cosine magnet to cancel uniform transverse mag- 
netic field components. These coils were also used to 
rotate the magnetic field for the transverse-polarization 
calibration runs. Additional trimming of non-uniform 
components of the magnetic field was necessary due to 
the magnetized steel in the floor and magnetic fields from 
other instruments in the NCNR guide hafl [SnillD|. The 
currents in all coils were nominally constant with no ac- 
tive compensation; however, all currents were monitored, 
and the fields in the vicinity of the decay region were 
continuously measured using two 3-axis fiux-gate mag- 
netometers. Magnet-current and detector-field monitor 
data were used to define cuts on the data as described 
section IVilAl 

The magnetic field in the detector region was mapped 
with a 3-axis fluxgate magnetometer before and after the 
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FIG. 5: Neutron polarization map showing the position de- 
pendence of the estimated lower-limit on neutron polarization 
assuming P — Ap. The top and bottom are cut off by the 
coUimating mask that defined the aperture of the analyzer. 
The increase of polarization from left to right is due to the 
PSM bending properties. 



emiT-II run. The large axial field of 560 f/T made pre- 
cise measurement of transverse component difficult due 
to flexing of the support rail for the magnetometer; how- 
ever, it was possible to estimate transverse components 
of 1-3 f/T. The corresponding magnetic-field misalign- 
ment was measured to be as large as 5.4 mrad, which 
was consistent with estimates based on the transverse- 



polarization calibration runs described in section VI E 



D. Detector overview 

The detector array is shown schematically in Fig. [2] 
The main detector chamber was milled from a single 
block of aluminum. The chamber was supported on kine- 
matic mounts and included mounts for cross-hair assem- 
blies at both the upstream and downstream ends. Beam- 
line and detector components were mechanically posi- 
tioned and aligned to better than 1 mm and 2 mrad, and 
the alignment was checked with crosshairs during final 
assembly. 

The vacuum during data taking was maintained us- 
ing a cryopump just downstream of the main detector 
chamber. In addition, liquid-nitrogen cooled cryopanels 
were situated at the ends of the proton detector assem- 
bly. The vacuum measured in this region was typically 
in the range of 3-4xl0~^ Torr. This low pressure signifi- 
cantly reduced backgrounds due to neutron interactions 
with residual gas and minimizes the possibility of scat- 
tering and neutralization of decay protons. On the other 



hand, changes in proton detector dead-layer of roughly 
3% per month, most likely due to the build-up of water 
and other volatile materials, were observed. 



E. Proton detectors 

The proton-detector cell is illustrated in Fig. [3] Each 
cell consisted of a grounded box with the top and the 
upper half of the sides covered by a grounded wire mesh 
(97% transmitting) through which the recoil protons en- 
ter. Once inside the box, the protons were accelerated 
and focused onto the SBD (Ortec model AB-020-300-300- 
S) by the field produced by a cylindrical tube maintained 
at a negative potential with respect to ground. During 
the course of the experiment, the acceleration voltage 
was varied in the range 25 kV to 31 kV. Geometric con- 
straints required that the SBD be positioned off-center 
relative to the focusing tube. In addition, the cryogenic 
epoxy holding the silicon wafer into its mount varied in 
thickness by up to one millimeter around the edge of the 
detector ring. These and other imperfections were in- 
cluded in the Monte-Carlo model of the proton focusing 
discussed in Section|Tn| The proton SBD active layer was 



300 /ini thick and 300 mm in area. 

The proton detectors were periodically calibrated in 
situ with ^'^^Am and ^°^Cd gamma-ray sources. Typical 
detector resolution was 4 keV FWHM but varied by sev- 
eral keV from channel to channel and over time. In most 
channels, thresholds were adjusted to reduce the count 
rates. This also resulted in truncating part of the ac- 
celerated proton spectrum introducing a proton-energy 
dependent efficiency and systematic effect discussed in 



section VI C Low energy tails on the proton peaks, ob- 
served in some channels, were of particular concern be- 
cause of the thresholds. The source of these tails was 
inconclusive and may have been due to a combination 
of a number of possibilities including dead-layer buildup 
and scattering from the focusing tubes themselves. In 
addition, reversible build-up of the detector dead layers 
was observed to occur over time. 

During the experiment detectors showed a variety of 
problems, including high leakage current and breakdown, 
and a portion of the data were taken with fewer than 
64 operating SBDs. Over the course of the experiment, 
however, only a few SBDs were not operational at any 
one time, and when averaged over the entire run, every 
proton cell had a duty factor greater than about 90%. 



F. Electron detectors 

Electron detectors are 50 cm by 8.4 cm by 0.64 cm 
thick BC408 plastic scintillators with Burle-8850 photo- 
multiplier tubes (PMTs) at each end. The sensitive re- 
gion of the electron detectors was fabricated from slabs 
of plastic scintillator cast to a thickness of 0.64 cm and 
diamond milled to a rectangular prism measuring 50 cm 
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by 8.4 cm. The thickness was sufficient to stop a 1-MeV 
electron and was therefore adequate for detecting neu- 
tron decay electrons, which have an endpoint energy of 
782 keV. The scintillators were wrapped in aluminized 
mylar and 20 fim thick aluminum foil to prevent charging 
and to shield the detectors from x-rays and field-emission 
electrons. This thickness of aluminum stopped electrons 
of energy up to about 50 keV at normal incidence and 
led to energy loss of 20 keV for incident electrons with 
energy 200 keV. The foil and mylar were included in all 
Monte Carlo simulations. 

Scintillation photons guided to either end of the scin- 
tillator by total internal reflection at the smooth surfaces 
were detected by the PMTs. An electron event requires 
coincidence of both PMTs within a timing window of 100 
ns. The time difference of the phototube signals at the 
two ends of each scintillator was available to determine 
the electron position with resolution of about 10 cm but 
was not used in the analysis. 

The electron detectors were calibrated periodically 
during the run using a ^"^Bi source. Conversion elec- 
trons from the source passed though a thin Kapton win- 
dow into the vacuum chamber and were incident on the 
backside of the scintillator slabs passing through a hole 
in the mylar and aluminum foils. Gain drifts of up to 
3% per month were observed, and PMT base voltages 
were occasionally adjusted to stabilize the gains. Re- 
maining drifts were compensated run-by-run in the anal- 
ysis using features of the neutron-decay-electron-energy 
spectrum. A threshold of SOkeV detected electron en- 
ergy was applied to the data in order to eliminate the 
effects of residual gain drifts. The effects of remaining 
small nonuniformitics in electron-energy threshold and 



response are discussed in Section VI C 



G. Monitors 

The following parameters were continuously logged 
during data taking: 

1. Magnetic fields at two positions within the detec- 
tor magnetic-field coils using 3-axis fluxgate mag- 
netometers 

2. Currents in all magnetic field coils including 
neutron-guide solenoids and longitudinal and trans- 
verse detector coils 

3. The neutron capture-flux monitor 

4. Proton-cell acceleration voltage 

5. SBD bias voltage and current. 

The cuts based on these monitor data were used to test 
for systematic effects and for cross checks of the results 
as discussed in section IVIIAl 
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FIG. 6: Two-dimensional spectra showing proton energy vs. 
proton time-of-flight vs. log of counts for all data. The boxes 
show the cuts used in the analysis (decay window) and to 
determine the random-coincidence background contribution 
(pre-prompt window). The pre-prompt events consist pri- 
marily of accelerated decay protons in coincidence with an 
uncorrelated electron-detector hit. 



V. DATA 

The neutron-decay data consist of proton-electron co- 
incidence events that record which detectors were hit, 
amplitude of the pulse in each detector, spin-flipper state, 
and proton time-of-flight (tep)- Data were acquired from 
October 2002 to November 2003 with some breaks for 
reactor and detector maintenance. The data were sepa- 
rated into runs of up to about four hours duration, and 
a total of 934 runs are included in the final data set. 

Fig. [6] shows the detected proton energy versus pro- 
ton time-of-flight for all events. Events prior to t^p — 
are random coincidences and are used to estimate back- 
ground rates as described in section |VI A| The fea- 
ture near t^p = primarily consists of true coinci- 
dences that are very closely spaced in time compared to 
neutron-decay-proton-electron events. These prompt co- 
incidences arise from several sources including cosmic-ray 
muons and neutron capture , which produces gamma rays 
that produce a false proton-electron coincidence through 
Compton scattering or pair production. In Fig. [7] we 
show the proton-ADC spectrum for the timing window 
—0.75 /xs < tep < 0.12 ^s as a function of detected SBD 
amplitude. Protons, primarily from neutron decay are 
accelerated in the focusing cells and produce the narrow 
peak at « 25 keV. The broad feature at higher energy is 
the minimum-ionizing peak for relativistic charged parti- 
cles, e.g. from high-energy gamma interactions and cos- 
mic rays. The pre-prompt data were used to estimate 
the random-coincidence contribution to the backgrounds 
discussed in section IVIAI 

Typical count rates were 3 s~^ and 100 s^^ for single 
proton and electron detectors, respectively, while the co- 
incidence rate for the entire array was typically 25 s~^. 
A total of 4.7 X 10® raw events were acquired. Of these, 
limits on magnetic field, leakage current, and other mon- 
itors removed 12% of the events, the analysis threshold 
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TABLE IL Summary of the final event selection. 



Description 



Approximate number 



Raw events 
Monitor cuts 

Electron-energy threshold 
Electron-detector multiplicity 
Final event sample 
Random coincidences 



4.7 X 10* 
5.6 X 10^ (12%) 
6.6 X 10^ (14%) 
3.3 X lO'^ (7%) 
3.2 X 10* 
1.0 X lO'^ 



on electron energy removed 14% and the single electron 
hit requirement removed 7% of the raw events. The data 
were filtered by cuts on the parameters described in sec- 



tion |IVG| as well as on the detected electron- energy and 
the requirement of a single-electron-detector hit. We also 
rejected events during a spin flip and during an unpaired 
spin-flip cycles at the end of a run in order to ensure 
equal time in each spin-flipper state. 

Most of the data were taken with nominal proton- 
acceleration voltages of 28 kV with smaller data sets at 
25 kV, 27 kV, and 31 kV. The 27-kV data were acquired 
during the initial running stages in 2002. Discussion of 
the data and results for different proton-acceleration volt- 
ages, running conditions, and detector configurations are 
discussed in section IVIIAl 



VI. SYSTEMATIC EFFECTS 



Eq. 



14 is based on an ideal symmetric experiment, 

a) accu- 



which incorporates the following assumptions: 
rate determination of the number of coincidences for each 
proton-cell-electron-detector combination, i.e. accurate 
background corrections; b) absence of proton and elec- 
tron backscattering; c) symmetry of the detector, specif- 
ically the equivalence of proton cells and uniform proton 
and electron detection efhciencies; d) uniformity of the 
neutron beam; e) uniform longitudinal polarization; and 
f) accurate determination of Ko and P. A number of 
effects break the symmetries of the experiment and are 
accounted for by corrections to D discussed in detail this 
section. 



A. Backgrounds 

A valid event is defined as the coincidence of a proton- 
detector signal and electron-detector signal that meet the 
selection criteria on proton and electron energy, proton 
time-of- flight, and the monitor-data. Electron detection 
requires a coincidence of both phototubes in one scintil- 
lator panel with detected electron energy above 80 keV. 
The proton-energy-time-of-flight window is shown in 
Fig. |6] Background events that satisfy these criteria arise 
from several sources. Estimates of the background frac- 
tions e^'*^^ for each component are discussed below, and 



Accidental Coincidences 

The count rates for each spin-flipper state were cor- 
rected for accidental coincidences on a run-by-run basis 
by scaling the counts in the pre-prompt timing window 
(-12.3 ^is < tep < -0.75 us) shown in Fig. [6] These 
random coincidences had the expected exponential de- 
pendence with respect to proton time-of- flight. A fit 
revealed a time constant of (2690 ± 730) /iS, which is 
sufficiently long that we made no adjustment for proton 
time-of-flight. This results in a correction on the back- 
ground subtraction estimated to be (2 ± 2) x 10^^. 



Prompt Coincidences 

An energy spectrum of events in the prompt window (- 
0.75 lis to 0.12 lis) is shown in Fig.jTj The dominant com- 
ponents are the peak at « 25 keV due to the accelerated 
protons from neutron decay, minimum ionizing charged 
particles predominantly due to cosmic-ray muons, and a 
background continuum. The low-energy behavior of the 
continuum was studied with high-voltage-off data, which 
eliminated the accelerated-proton peak. This component 
appears to fit well to a double exponential. 



TABLE III: Averages over all detectors of background frac- 
tions for each background component. 



Source 


It (10-^) 


Accidental coincidences 


2±2 


Low-energy continuum 


6± 2 


Minimum Ionizing 


0.003 ±0.002 



Background-related contributions to D 

The asymmetries w^^'^^ must be corrected to determine 
the background- free asymmetries lUg'*^^. The correction 
is 



PiSj PiEj 



(15) 



the averages over all detectors, e;,, are given in table III 



The first term is the multiplicative correction to lyP'^j 
due to the dilution of the asymmetries by backgrounds. 
This dilution can produce a contribution if the back- 
ground corrections are non-uniform and do not com- 
pletely cancel when the w^^'^^ are combined into D. The 
multiplicative correction to D was determined using mea- 
sured e^'^^J and w'^^'^^ according to Eq. 15 resulting in a 
correction to b of (0.03±0.09) x 10"''. The multiplicative 
correction is dominated by the exponential components 
of the prompt events shown in Fig. [7| The second term 
in Eq. 15 is due to a possible asymmetry in the back- 
ground, w^''^^. The w^'^^ were found to be uniform and 
consistent with zero ac ros s the detector and were com- 
bined according to Eq. 14 The resulting Db was scaled 
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FIG. 7: Typical proton-ADC spectrum for the "prompt" tim- 
ing window —0.75 ps < t^p < 0.12 jis. The lines show best-fit 
contributions due to neutron-decay protons, minimum ioniz- 
ing, and an exponential tail determined from runs with zero 
proton-acceleration voltage. 



by the sum of for all contributions to determine the 
additive correction to D of (-0.07±0.07) x 10"^. For the 
result reported in reference [TT] , no correction was made 
for the additive-background effect. 



B. Backscattering 

Backscattering events fall into two categories: 1) par- 
ticles scattered from somewhere in the apparatus to the 
proton or electron detector, and 2) incident protons and 
electrons that scatter from the respective detectors with- 
out registering a hit. Particles backscattered from the 
proton and electron detectors lead to multiplicative and 
additive corrections to the w^' that can affect D in a man- 
ner similar to background effects. B acksc attering also 
affects and is discussed in section 



VI F 



Backscattering corrections to D were determined by 
Monte Carlo simulations with empirical validation based 
on studies of the 35° and 55° proton-electron coinci- 
dences. Backscattering probabilities are similar for all 
proton-clectron-detector combinations; however, the rate 
of neutron-decay coincidence counts for the 35° and 55° 
pairings is about a factor of ten less than for the 125° 
and 145° pairings used to determine D. Data from the 
35° and 55° pairings are shown in the electron-energy- 
proton-time-of-flight plane in Fig. [Sj The boxes indicate 
three distinct kinematic regions: primary neutron-decay 
events for the 35° and 55° proton-electron coincidences 
are restricted to a region included in region A; region B 
consists of higher energy electrons, which are primarily 
electrons with initial momentum directed at large an- 
gles with respect to the detected proton but detected 
in the 35° or 55° electron detector due to backscatter- 



FIG. 8: Data for electron energy vs. proton-electron delay- 
time for proton cell pi in coincidence with electron detector 
ei at 55° . The boxes indicate regions of kinematically allowed 
proton-electron coincidences (region A) , electron backscatter- 
ing from detector elements across from the proton cell (region 
B), and proton backscattering (region C). 



ing; region C corresponds to high energy electrons from 
events with large proton-electron angular separation and 
delayed proton time-of- flight, which are primarily due to 
proton backscattering. The backscatter fractions are es- 
timated from the data by integrating within the boxes 
in Fig. [8] or by fitting the time-of-flight spectra and are 
consistent with the Monte Carlo predictions. 



Electron backscattering 

Electron backscattering was incorporated into the 
Monte Carlo simulations using penelope. The multi- 
plicative correction to D was determined using estimates 
of the backscattering fractions from the Monte Carlo to 
correct the measured w^'^^^ similarly to the background 
corrections given in Eq. |15| With the correction ap- 
plied, D changes by (0.11 ± 0.03) x 10~^, where the un- 
certainty reflects the 20% uncertainty assigned to the 
Monte-Carlo results due to limitations of the detector 
and beam model and limited knowledge of backscattering 
at these low energies. The additive correction would van- 
ish for isotropic backscattering and uniform detector effi- 
ciency; however, the beam expansion, magnetic field, and 
detector elements break these symmetries. The result for 
the additive correction to D is (0.09 ± 0.07) x lO^"', and 
is limited by the statistical precision of the Monte Carlo 
due to the small fraction of backscatter events. This 
additive-electron-backscattering correction was not ap- 
plied in the analysis reported in reference |11| . 
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Proton hacks cattering 

Backscattered protons in principle affect the asym- 
metries similarly to backscattered electrons. Proton 
backscattering probabilities determined with SRIM were 
input into the Monte Carlo. While the proton backscat- 
tering probabilities at the energies characteristic of neu- 
tron decay are comparable to electron backscattering 
probabilities, there is a large probability for neutral- 
ization [35]. Detailed calculations lead to the conclu- 
sion that proton-backscattering effects on D are approx- 
imately an order of magnitude smaller than electron- 
backscattering effects, and we set an upper limit on the 
correction to D of less than 0.03 x lO""*. This upper limit 



is reflected in the uncertainty given in Table IV 
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C. Efficiency corrections 



FIG. 9: Results for v^"" averaged for each of the four proton- 
detector planes. The dashed line is a fit to sin(/)^», and the 
solid line includes a sin20''' term. Monte-Carlo simulations 
show that a sin 20''' term will arise from « 10 keV differ- 
ences of the electron-detection thresholds for upstream and 
downstream ends of the electron detectors. 



The wP'^^J from Eq. [9] depend directly on the energy 
dependence of the efficiencies through the beta and neu- 
trino correlations. Spatial variations of the proton and 
electron energy dependence break the symmetry assumed 
in combining proton-cell data as given in Eq. [14] In this 
section, we discuss the effects of axial and azimuthal vari- 
ation of the electron energy thresholds and the variation 
of the proton-energy-dependent efficiency as a result of 
changes of individual SBD thresholds. 



1. Nonuniform electron- energy thresholds 

From Eq. [12] it can be seen that the proton-electron 
correlation leads to a difference of /3e for the two electron 
detectors that pair with pi {e.g 63 and 62). This leads 
to a contribution to u^* from the bet a- asymmetry {A- 
coefficient correlation), which cancels for a uniform neu- 
tron beam and detector when the u^' from two axially- 
symmetric proton cells are averaged {e.g. cells pi and 
P15). In the case of nonuniform electron-detector efficien- 
cies, specifically due to spatial variation of the electron- 
energy threshold, this cancellation is not perfect and 
leads to a dependence on the azimuthal proton-cell po- 
sition, i.e. a dependence on proton-detector plane; how- 
ever, another symmetry arises because each electron de- 
tector is paired with two azimuthally opposed proton 
cells and contributes to the two w^' with opposite sign. 
This results in cancellation of any efficiency variation to 
first-order in siiKpP', as given in Eq. 13 The sin20^'' 
dependence of however, does not cancel. 

The electron energy thresholds were measured to dif- 
fer between upstream and downstream phototubes by as 
much as 10-20 keV resulting in the sin2(/)P' dependence 
evident in Fig. [9] where the combined data for w^*' for all 
proton cells in a single detector plane are shown as a func- 
tion of . The dashed line shows the best fit to sin (j>P^ 
only, and the solid curve includes a sin 20^' contribution. 



Monte-Carlo studies with an upstream-downstream vari- 
ation of electron threshold of approximately 10 keV for 
a single electron-detector show a similar sin 2(/)P' contri- 
bution. 

The Monte Carlo also shows that the axial depen- 
dence of the electron thresholds are cancelled when all 
four proton-detector planes are combined and therefore, 
does not contribute significantly to D as long as the pro- 
ton efRciencies are energy independent. The variation 
of proton energy dependent efficiency, discussed below, 
could couple with an electron-threshold variation. Cou- 
pling between the observed relative efficiencies of proton 
detection and the Monte-Carlo model of electron thresh- 
old variations had an effect on D of (— 0.04±0.1) x lO^''. 



2. Proton-SBD efficiency variations (Shift-threshold effect) 

A set of typical calibrated proton-cell SBD spectra is 
shown in Fig. [lO] for both spin-flipper states and for both 
125° and 145° coincidence angles. The electron detector 
at 145° has approximately twice the rate as the 125° de- 
tector due to the proton-electron correlation, which fa- 
vors back-to-back proton and electron momenta. The 
spectra for the two spin-flipper states differ in both area 
(total number of counts, N^*^' ) and the spectral distribu- 
tion due to the beta asymmetry, which affects the angular 
distribution of electrons and thus the proton momentum. 
The energy deposited in the SBD in turn depends on the 
incident proton momentum, position and angle of inci- 
dence on the proton cell. We characterize each spectrum 
by a centroid, which depends on the neutron spin sate. 
The difference of the centroids for the two spin-flipper 
states varies with the axial position of each proton cell 
and varies in magnitude from zero to about 125 eV. In 
Fig. [11] we show the difference of the centroids as a func- 
tion of proton-cell position for a single proton-detector 
plane (pi-pie) paired with the electron detectors at 125° 
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FIG. 10: Proton-cell SBD spectra for a single proton detector 
(e.g. pi) for both electron detectors 62 and 63 for both spin- 
flipper states: up (neutron polarization parallel to B) and 
down (neutron polarization anti-parallel to B). 



FIG. 11: Differences of proton-detector centroids for proton 
cells P1-P16 as a function of axial position for both 125° and 
145° proton-electron-detector pairs. Fluctuations from cell- 
to-cell reflect variations of energy loss and other features of 
individual proton cells. Odd and even proton detector refer 
to pi, Pa, etc. and p2, P4, etc., respectively. 



and 145°. 

For the SBD spectra shown in Fig. [TOj the low energy 
portion of the spectrum shows that the applied thresh- 
old allows detection of almost all of the protons for both 
electron detectors and both spin-flipper states. The SBD 
thresholds were adjusted throughout the run in order to 
maintain the electronics' noise at a manageable level and 
minimize dead-time losses. In some cases, the thresh- 
old cut significantly into the SBD spectra introducing 
a proton-energy-dependent detection efficiency. Due to 
the dependence on neutron polarization, this energy de- 
pendence could lead to a significant effect on the wP-'^j 
for an SBD with a high threshold, though the effect on 
tif' is largely mitigated because the low-energy portion 
of the proton energy spectrum is nearly equally affected 
for the two electron detectors. Fig. [12] shows that there 
are several notable anomalies of the for example 

for proton cells p4 and P12; however, as expected, the 
anomalies are similar for both electron-detector pairings. 
Thus the , as shown in Fig. [13) do not appear to fluc- 
tuate significantly compared to other proton cells. 

To estimate the effect on 13, neutron-polarization de- 
pendent proton energy spectra were generated by Monte 
Carlo for all proton-detector-electron-detector pairs and 
convoluted with model detector response. The model de- 
tector response (modified gaussians with separate widths 
above and below the centroid) were based on measured 
proton-SBD spectra. These spectra and the thresholds 
varied significantly during the course of the experiment, 
and the set of response functions was therefore based on 
an average over the data subsets. Parameters of the re- 
sponse functions were varied over a range characteristic 
of the variations during the experiment in order to esti- 
mate the uncertainties. The estimated correction to D is 
(-0.29 ± 0.41) X 10""^. The effect was also estimated by 
correcting uf^'^^ on a day- by-day basis using the centroid 
shift and slope at threshold with consistent results. 
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FIG. 12: ■w^^'^^ for all proton cells paired with electron de- 
tectors at ±125° and 145° with respect to the proton plane. 
The open triangles pointing up are for positive cross product 
Pp X Pe, and the filled triangles pointing down are for negative 
cross product. The u^* shown in Fig. [13] are found from the 
difference of the two w'''^'^^ for each proton cell. 



D. Beam Expansion Correction 

The neutron beam expansion affects the cancellation 
of the beta-asymmetry (Eq. 12 ) for axially-paired pro- 
ton cells because the average of the electron velocity (/3e) 
depends on the radial size of the beam: for decay ver- 
tices further from the center of the detector, the proton- 
electron angular separation is larger, which corresponds 
to higher energy electrons. This would be largely can- 
celled by the combination of the v^' prescribed by Eq. [llj 
because the difference of (/3e) enters with opposite sign 
for two adjacent proton cells (pi and P2); however, the 
magnetic field also affects the electron-proton angular 
correlation and therefore (/3e) resulting in a magnetic- 
field dependence. The 560 /xT magnetic field makes the 
magnitude of the difference for the two electron detec- 
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FIG. 13: 11^' for all proton cells for the experiment. The 
solid and dashed lines are from the Monte Carlo simulation 
for realistic experimental conditions shown in Fig. [l4]with ad- 
justment for each proton-cell plane to account for transverse 
polarization and beam-shape effects. 



tors {e.g. pie^ and ^162) smaller for odd detectors (pi, 
P3, p5, etc.) than for even detectors {p2, pa, pe, etc.). 
The result is that the slopes of u^* vs. z^' have different 



magnitudes for even and odd proton cells. Fig. 14 shows 
a comparison of a uniform, nonexpanding neutron beam 
and the realistic beam based on the beam-distribution 
maps of Fig. |4j A possible systematic effect arising from 
the beam expansion coupling with a polarization gradient 
is discussed in section IVl El 

To estimate the correction, Monte Carlo simulations 
were run using the measured neutron-beam distribution. 
A number of simulations were run with different beam 
distributions produced by shifting the maps by up to 2 
mm and rotating by 5°, both significantly greater than 
the mechanical alignment constraints on the beam-line 
and detector components. In Fig. [15] we show the results 
for the beam expansion effect as a function of magnetic 
field. The estimated correction at _B = 560 /iT is D 
is (—1.5 ± 0.4) X 10"''', where the uncertainty is mainly 
due to the change of the effect as the beam-distribution 
models were varied. 



E. Transverse polarization effects 



As discussed in section 11 B transverse asymmetries 



arise in the experiment due to the polarization misalign- 
ment with respect to the detector axis. The misalignment 
was determined from the data combined with transverse- 
polarization calibration runs combined with Monte-Carlo 
simulation of the beam-shape effects and by magnetic 
field maps measured before and after the run. The trans- 



u!j^' given m 



verse polarization leads to a contribution 
Eq. 13 For uniform-symmetric beam and polarization, 
cx sin(0p — and is cancelled when data from 
two opposing proton cells are combined; however, due 
to the asymmetric emiT-11 beam, the cancellation is not 
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FIG. 14: The broken symmetries due to magnetic field and 
beam expansion contribute to a D for longitudinal polariza- 
tion. The 560 pT magnetic field makes the magnitude of u*"' 
larger for even detectors than for odd detectors, i.e. the de- 
pendence of v^^ vs. z^^ differs largely due to the difference 
of for the two electron detectors, e.g. the expanding 

beam makes the magnitude of ((/3e^''^) — {Z??^*^^)) larger for 
downstream detectors since the larger beam leads to larger 
proton-electron angular separation. 
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FIG. 15: The beam expansion effect as a function of mag- 
netic field determined by Monte-Carlo simulations using the 
measured beam profiles. The solid line is a fit to the slope, 
constraining the effect to zero at 5 = 0. The error bars rep- 
resent the statistical uncertainty of the Monte-Carlo. 



complete. Displacement of the beam perpendicular to 
the transverse polarization would result in a contribution 
to D that is the product of the polarization misalign- 
ment and the perpendicular misalignment of the beam. 
This is called the ATP or asymmetric-beam-transverse- 
polarization effect. An additional effect arises if the 
transverse polarization is nonuniform, specifically that 
sin (/)p varies along the axis of the detector. 



1. Asymmetric Beam/Transverse Polarization Effect 

We have studied the ATP effect using transverse po- 
larization calibration runs. Data from a transverse- 
polarization calibration run with {0p,4>p) = (90°, 180°) 



is shown in Fig. 16 This amplifies the ATP effect, which 
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FIG. 16: Results of a transverse-polarization calibration run 
with 4>p = 0°, I.e. polarization pointing in the x direction 
(see Fig. [2] for coordinate system definitions). Each point 
combines data for two proton cells with the same \z°'\. The 
solid line is a fit to the form of Eq. |13| The average amplitude 
for 11 calibration runs is < vt >= 0.4565 ± 0.0128, where the 
uncertainty is the standard deviation of the 11 amplitudes. 



is then scaled by the size of the transverse polarization 
for the run to determine the effect on D. Calibration 
runs at several azimuthal angles (j)p map out the effect 
to provide a more accurate estimate. The size of the ef- 
fect was also estimated in Monte Carlo simulations with 
Op — 90° and several different values of (j)p. 

For calibration runs, the axial magnetic-field coils were 
used to cancel the Earth's field component along the de- 
tector axis, and the transverse-field coils were used to 
produce a field of about 100 /iT perpendicular to the 
detector axis. The azimuthal angle could be selected; 
however, due to power supply limitations, the magnitude 
of the transverse field could not be maintained at 100 /iT 
for all azimuthal angles. In Fig. 16 we show the average 



vfj: for even and odd proton cells as a function of from 
a transverse polarization calibration run with Op — 90° 
and (/)p = 0°. A total of eleven polarization-calibration 
runs were taken over the course of the experiment. For 
each run, the amplitude and the offset are determined by 
fitting the data to a sinusoid. The average amplitude for 
all calibration runs with 6'p = 90° is 



< Vt >= 0.456 ±0.013, 



(16) 



where the uncertainty is the standard deviation of the 
eleven amplitudes. 

For the data runs, Op is small; thus the asymmetric 
beam shape contributes significantly to the (/)^' depen- 
dence of the vs. The beam maps shown in Fig. [4] show 
that, in addition to the beam expansion, the center of 
gravity of the beam is horizontally displaced from the de- 
tector axis and rises in the vertical direction for the down- 
stream map. To separate the beam-shape and transverse- 
polarization effects, the w^' were corrected for the beam- 
shape effect, which was determined from Monte Carlo 
simulations with no transverse polarization. The value 
of sinOp is found by dividing the amplitude of a sinusoid 
fit to the corrected u^' by the average amplitude from the 



transverse-polarization calibration runs. In Fig. 17 we 
show results of Monte-Carlo studies for the beam-shape 
along with data and the resulting estimate of the trans- 
verse polarization effect. The Monte Carlo results for the 
beam-shape effect are very sensitive to the small differ- 
ences in the beam-map registration and orientation. This 
results in large uncertainties in the transverse polariza- 
tion direction. The results are 

sin Op = (8.5 ± 4.3) x 10"^ (j)p = (40 ± 72)° . 

The uncertainty is determined from the uncertainty on 
{vt) and errors arising from the beam-shape correction 
and the sin2(/)p effect (Fig. |9|. Alternatively, using the 
magnetic-field maps discussed in section |IVC| we esti- 
mate 

sin 6lp = (3.8 ± 1.0) X 10"^ (/)p = (45 ± 5)°. 

which is consistent with the analysis using the data 
and beam-shape correction. To determine the ATP 
correction, we use the 1-a limits from the transverse- 
polarization calibration runs, which give 

Op < 12.8 X 10"^ - 32° <<j)p < 112. 

These limits provide more conservative bounds on the 
ATP effect than those base on the magnetic-field maps. 

we show the data with fit and the re- 



in Fig. 18 



suit for the Monte-Carlo model of the experiment for 
sin 6'p = 8.5 x 10~^. From the fit and the values of Op 
and tpp given above, we determine the correction to D of 
(-0.07 ±0.72) X 10-"*. The relatively large uncertainty 
is due to the uncertain azimuthal orientation of the po- 
larization. 



2. Polarization twist 

A twist of the polarization, that is an upstream- 
downstream difference in 0p, would not be cancelled in 
combining 16 proton cell w^' into D. A limit on the 
polarization twist was estimated by separating data for 
upstream (z^' < —6 cm) and downstream (z^"' > 6 cm) 
proton cells and correcting for beam-shape effects. The 
beam-shape corrected data are shown shown in Fig. |19| 
Sinusoid fits provide estimates for (pp for upstream and 



downstream portions of the detector of < 



14.3°±5.2° 



and 



idn 



= 16.6° ± 6.0°, from which we use to estimate 



a maximum upstream-downstream difference A(f>p = 
13.5°. We double A(j>p to account for a change over the 
entire length of the proton-detector plane and use the fits 
from Fig. [18] to set an upper limit on the magnitude of 
the ATP-twist effect of D of 0.24 x 10"''. This upper 



limit is reflected in the uncertainty listed in Table IV 



3. Polarization nonuniformity 

A nonuniform neutron polarization affects the cancel- 
lation of the electron (A) and antineutrino (B) asymme- 
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FIG. 17: Beam-shape and transverse-polarization effects on 
v^^ . The filled squares show the average v^^ for each proton- 
detector plane as a function of azimuthal angle for the exper- 
iment. The solid line is a fit to the data including the sin 20^* 
contribution. The dotted and dashed lines are, respectively, 
the estimated effects of the beam shape and the transverse po- 
larization. The beam shape effect was determined by Monte 
Carlo simulations for the beam maps shown in Fig. |4] the 
transverse polarization was estimated from the difference of 
the sin (^'''-only fit (Fig. [9| and the beam-shape effect. 
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FIG. 18: ATP effect on D as a function of <f>p . The data points 
were determined from the eleven polarization-calibration runs 
scaled by sinOp = 8.5 X 10"^. The solid line is a sinusoidal fit 
to the data with the offset constrained to be zero. The dashed 
line, given by 0.48 x 10"" sin((?!)p - 204°), is the Monte-Carlo 
prediction based on beam maps and sinOp — 8.5 x 10~^. 
The horizontal bar extending from -33.1° to 112.1° shows the 
range of (jjp based on Monte Carlo estimates using reasonable 
variations of the beam-shape map registration. 



tries in D. The polarization map, which measures the 
combination of neutron-beam polarization P and the an- 
alyzer power Ap, is shown in Fig. [5] Assuming P = Ap, 
the polarization as a function of position was convolved 
with the density of neutron-decay vertices determined 
by Monte Carlo to determine an effective polarization 
for the decays detected by each proton cell. This was 
also studied using data from both transverse-polarization 
calibration runs and normal running. This effective po- 
larization varied by a few percent along the beam axis 
due to the beam expansion; however, the average over all 



FIG. 19: Data, corrected for beam-shape effects, for v^"- for 
upstream and downstream proton cells. The solid lines in- 
dicate sinusoidal fits with (f)^ = 14.3° ± 5.2° and 
16.6° ±6.0°. 
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cells in each proton-detector varied by less than 0.005. 
The effect on D was estimated by Monte-Carlo using the 
averaged polarization for each proton-detector plane. A 



possible transverse polarization gradient or ^) was 
also investigated with the Monte Carlo and found to have 
a negligible effect. 



Uncertainties in Polarization and Ko 



Uncertainties in P and Kpi lead to errors proport ional 



to D. The polarization analysis described in section IV B 
results in the lower limit P > 0.91 (90% c.L), and we take 
P = 0.95±0.05 for the purpose of analysis. The resulting 
contribution t o th e uncertainty on D is 0.05D, which is 
given in Table IV The kinematic quantities depend 



on the proton- electron angular separation due primar 
ily to the cross-product Pp x Pe, and are different for 
the 125° and 145° proton-electron-detector pairs. The 
K'^ were determined by Monte Carlo simulations, and 
a number of parameters were varied including beam size 
and shape, electron backscattering fraction, and electron- 
energy threshold. For proton-detector-electron-detector 
pairings at 125° and 145° we find \K}^^ \ = 0.420 ±0.008 
and l^^jj^^l = 0.335 ±0.016, where the uncertainties esti- 
mate the variation of Kp) as the Monte Carlo parameters 
are varied. The average, used to determine D (see eq. 12 ) 
is 



Kn = 



\K 



D 



= 0.378 ±0.019. 



(17) 



The resulting contribution to the uncertainty is 
0.051) = 0.04 X 10-*. 



G. Spin flip and timing errors 

Errors in the w^^^^ can also arise if the neutron polar- 
ization, flux, or counting time in each spin-flipper state is 
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dependent on spin-flipper state. In these cases, the error 
is approximately proportional to D. For example if the 
neutron flux or density within the detector is dependent 
on spin-flipper state, the count rates would be modified 
from the such that N^^' = N^^' (1± %^), where N^'''' is 
the count rate for a spin-flip symmetric experiment and 
Sp is the fractional change in the flux. In the event of 
a dependence of the flux on the spin-flipper state, the 
asymmetries would be modified from the spin-flip sym- 
metric asymmetries Wq'''\ The measured asymmetry is 

^P.e, ^ ^P.e, ^p,e,^2^ ^ Sf{wI''^)\ (18) 

When the differences of the wP^"^^ are combined into the 
yP' and finally D, the correction is 

Do-D^SlD~0{D^), (19) 

where we assume the measured D is approximately equal 
to the corrected Dq. A spin- flip correlated counting time 
difference (St) and polarization difference (Sp) introduce 
similar corrections. 

Measured limits on the flux and counting time cor- 
relations with spin-flipper state are Sp < 0.004, and 
6t < 10~* respectively. Thus the corrections to D due 
to spin-flip correlated flux and counting time variations 
are less than 0.01 x 10^"*. The spin- flipper efficiency was 
estimated to be greater than 0.95, implying Sp < 0.05 
and a correction to D less than 0.04 x 10""*. 



TABLE IV: Systematic error corrections in units of 10^''. We 
determined upper limits on the magnitude of corrections for 
proton backscattering, polarization non-uniformity, and the 
ATP-twist effect, and thus these corrections are indicated as 
with the upper limit indicated by the uncertainty. Corrections 
for spin-correlated flux and spin-correlated polarization are 
less than 0.01 x 10~^, thus no correction was made, and the 
contribution to the uncertainty was negligible. 
^ Polarization and Kd are included in the definition of D. 
Assumes polarization uncertainty of 0.05. 

Source Correction Uncertainty 



* In reference |11| this entry had a typographical error. 



Background additive 


-0.07 


0.07 


Multiplicative* 


0.03 


0.09 


Electron backscattering additive 


0.09 


0.07 


Multiplicative 


0.11 


0.03 


Proton backscattering 





0.03 


Electron threshold non-uniformity 


0.04 


0.10 


Proton-threshold effect 


-0.29 


0.41 


Beam expansion 


-1.50 


0.40 


Pol. non- uniformity 





0.10 


ATP - misalignment 


-0.07 


0.72 


ATP - Twist 





0.24 


Spin-correlated flux 





<0.01 


Spin-correlated polarization 





<0.01 


Polarization 


a 


O.Ol"" 


Kd 


a 


0.04 


Total corrections 


-1.66 


0.97 



H. Summary of Systematic Error Corrections 



The corrections to D are summarized in Table lIVl The 
polarization and instrumental constant P and Kp) are 
included in the definition of D and are not included as 
corrections; however, the uncertainties on both are in- 
cluded in the table. The total correction is the sum of all 
corrections. The uncertainties are independent and are 
therefore combined in quadrature to determine the total 
systematic- error uncertainty. 



VII. RESULTS 



parameters with windows set by typical variations that 
produced useful data {e.g. for the flux-gate magnetome- 
ter z component: 5.56 G < < 5.62 G). To inves- 
tigate the effects of the cuts, the nominal window was 
expanded by a factor of two in most cases {i.e. 5.53 G 
< Bz < 5.65 G), and the change in D was noted. For the 
electron-energy threshold the cut was lowered to 40 keV, 
which was below the threshold in some channels, and for 
electron-detector multiplicity, the cut was changed from 
one to two. With the exception of the proton acceler- 
ation voltage discussed below, D changed by less than 
0.1 X 10^^ or 5% of the statistical error. 



A. Cross Checks 

Several cross checks have been performed to validate 
the analysis and search for systematic errors. All cross 
checks were performed on blinded data. The cross checks 
fall into three main categories: 1) varying the cuts on 
experimental parameters, 2) breaking up the data into 
subsets, and 3) alternative definitions of D that would 
be equivalent for an ideal experiment. 



Data Cuts 



The original cuts on experimental parameters listed in 



section IV G were established using nominal operating 



2. Data-subset studies 

Each independent measurement of D combines v^' 
from sixteen proton cells, all of which have relatively 
high efficiency for proton-electron coincidences; however, 
during the experiment, individual proton-cell SBDs were 
turned off for extended periods due to high leakage cur- 
rents or noise. Possible variations of the results due to 
varying experimental conditions were studied by break- 
ing up the experiment into subsets with roughly equal 
statistical weight. Subsets were separated by several pos- 
sible changes including proton-acceleration voltage, num- 
ber of live SBDs, and changes to the magnetic field prior 
to transverse-polarization calibration runs. For each of 
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FIG. 20: Results for D by data subset. The proton- 
acceleration voltage and the proton-cell sets included for each 
subset are indicated. (Proton-cell set 1 refers to \zp^ — ±2| 
cm, etc.) The weighted average of all subsets is 0.58 ± 2.14 
with = 10.44 for 12 de grees of freedom. 



the four I?, 16 operating SBDs are required; however, a 
number of data subsets had one or more sets of proton- 
cells missing, in which case we report the weighted av- 
erage of all available full sets of 16 proton-cells. The 
results for the subsets are shown in Fig. |20j A possible 
correlation of D with high voltage was revealed in the 
study of cu ts and investigated with data subsets shown 



in Fig. 
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Assuming D is independent of acceleration 
voltage results in = 10.4 for 12 degrees of freedom. 
Allowing a linear dependence of D with acceleration volt- 
age results in = 5.6 for 11 degrees of freedom. The 
change in implies a 2.1-sigma slope. In addition, the 
acceleration-voltage dependence of the focusing proper- 
ties was extensively studied by Monte Carlo simulations, 
which showed no effect correlating D with acceleration 
voltage. We therefore consider the preference for a slope 
to be an accidental correlation. 



3. Alternative definitions of D 
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averaged the small- 
which cancelled the 



The definition of D given in Eq 
est number of proton cells {i.e. 16) 
beta and neutrino asymmetries in the presence of trans- 
verse polarization and the 560 /iT magnetic field. Our 
final result, the weighted average of the four indepen- 
dent determinations of D, thus has the smallest possible 
statistical error. For a uniform detector and beam, D 
can also be defined as a) the simple average of all 64 v'^^ , 
which differs from D only due to small changes in weight- 
ing of the individual v^* ; b) the average of the combined 
data from four individual proton-detector planes; c) the 
offset when the w^' are fit to a sinusoid (Eq. 131; d) the 
offset for the sinusoid fit when the averages of f^' for 
each proton-detector plane; or e) the paired proton-cell 
approach used in the analysis of emiT-I data [9j. The 
results for D are all found to be consistent with D for 
all analyses based on w^', i.e. a)-d). The paired-proton- 
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FIG. 21: Results for D for the entire experiment for each set 
of 16 proton cells. The weighted average is 0.72 ± 1.89 with 



X 



0.8 for 3 degrees of freedom. 



cell approach, method e), is known to be very sensitive 
to the proton-threshold effect discussed in section |VI C| 
and produced a significantly different value. Correcting 
for the estimated proton-threshold effect on D based on 
day-to-day correction of the w^''^^ yielded a value consis- 
tent with D. Cross checks were also performed on data 
subsets and found to be consistent with the exception of 
the 31-kV data for which a large number of proton cells 
were dead. 

A blind analysis of the asymmetries was adopted by 
adding a quantity K^'^'B to each iu^-'^j (equation |9]) so 
that when D was extracted from Eq, 
from the true value by S, where —0.01 < B < 0.01. 



it was offset 
The 

factor B was revealed and subtracted as the final analysis 
step, after the corrections for systematic errors and all 
uncertainties were determined. 



B. Final Result 

When averaged over the entire run, each SBD was live 
for a majority of the time and had a high average effi- 
ciency. We can therefore combine counts for the entire 
run to determine the u^'s and to extract D for each set 
of 16 proton cells. The results for the four separate D 
are presented in Fig. [21] and Table |V] The weighted av- 
erage of the four D is (D) = (0.72 ± 1.89) x 10"''. When 
combined with the total of corrections listed in Table [TV] 
our final result is 



D = [-0.94 ± 1.89(stat) ± 0.97(sys)] x 10" 



(20) 



This differs from reference '11' due to refinement of the 
additive corrections due to background and backscatter- 



ing, which changed by 
spectively. 



-0.07 X 10-4 and 0.09 x lO"'^, re- 



VIII. CONCLUSIONS 

Our result represents the most sensitive measurement 
of D in nuclear beta decay and can be interpreted in 
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FIG. 22: Sensitivity to T violating couplings due to scalar and tensor currents from this work compared to limits from other 
experiments. Left: sensitivity to the iinaginary part of scalar couplings, assuming tensor currents to be purely real and equal 
to the largest presently-allowed value = 0.12. The bounds labeled "R from n" are from Ref. [32j, and those labeled "little 
a in 0^ — > 0^" are from a combination of Ref. [15] and [IQ. Right: sensitivity to the imaginary part of tensor couplings, 
assuming scalar currents to be purely real and equal to the largest presently-allowed value Cg = 0.1. The bounds labeled 
"little a from ®He" are from reference [47] and those labeled "R from *Li" are from Ref. ,484. The current work is labeled 
"emiT." All allowed regions represent 95% confidence levels. 



TABLE V: Results for the experiment for four proton-cell sets 
in units of 10~*. Uncertainties are statistical errors only. 



Proton-cell set 
1 
2 
3 
4 

Average 



\z\ (cm) 
2 
6 

10 
14 



D 

-2.57 ±3.49 
-1.57 ±3.67 
1.60 ±4.08 
0.20 ±3.93 
0.72 ± 1.89 



terms of possible extensions of the Standard Model. 
Rewriting Eq. [7] 



2|A| 



-3|AP 
sin (l)Av 



-Im{c'+(C+)*+Cs((7y)*} 



(21) 



where = {Cs ± C^)/Cy and = (Ct ± C'j,)/Ca, 
and we have neglected the order-a radiative corrections, 
which yield negligible contributions given the existing 
limits on scalar and tensor currents. 

Assuming no scalar or tensor currents, our result 
constrains the complex phase between the axial-vector 
and vector currents {Ca/Cv = \X\e^'^^^ ) to (j)AV — 
180.012° ± 0.028° (68% confidence level). If all cur- 
rents are allowed, for example due tor leptoquark ex- 



change, the equation contains five phases, making it dif- 
ficult to compare the sensitivity of our experiment with 
respect to other probes without further assumptions. In 
the specific case where there are no special cancellations 
between terms, we estimate the sensitivity of our mea- 
surement under two different assumptions. Fig. |22] left 
panel, shows the limits on the imaginary component of 
the scalar currents, assuming that the tensor currents 
are purely real and equal to the largest value allowed by 
present constraints. In this case, because the limits on 
are much smaller than those on our result is most 
sensitive to the imaginary component of Cg . The right 
panel shows a similar plot but now assuming that scalar 
currents are purely real and equal to the largest value 
allowed by present constraints. 



A. Potential Improvements 

The result for D presented here has comparable sta- 
tistical and systematic uncertainties. Thus an improved 
experiment with the same apparatus would need both 
more neutron decays and reduced systematic effects. A 
new beam-line (NGC) under construction at NCNR and 
the PF-1 beam at ILL could provide a factor of 10 or more 
decays per unit time. Reducing the three major system- 
atic corrections requires eliminating the proton-threshold 
variations, a more symmetric neutron beam, and smaller 
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magnetic field. The symmetry of the neutron-beam is 
most strongly affected by the supermirror-bender neu- 
tron polarizer, while the 560 //T magnetic field was cho- 
sen to effect sufficient velocity averaging of transverse- 
neutron polarization produced in the current-sheet spin 
flipper. An alternative polarizer is a steady-state po- 
larized ■^He spin filter [15]. Intense cold neutron beams 
have been shown to affect the rubidium and '^He po- 
larization |50j : however, this appears to be a solvable 
technical challenge [51]. The 560 fiT guide field can be 
reduced by using an adiabatic-fast-passage neutron spin 
fiipper and effective shimming of the magnetic field along 
with shielding of external field perturbations. Thus a 
factor of three or more improvement in sensitivity to D 
appears within reach with the current apparatus. Ex- 
tending the sensitivity to the level of final-state-effects 
[Dpsi ~ 10^^) and beyond is a well motivated goal that 
would require an apparatus with greater geometric effi- 
ciency for both proton and electron detectors. 
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